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Abstract. Modified Newtonian Dynamics (MOND) is evolving from an empirical to a decent theory respect- 
ing fundamental physics after Bekenstein (2004) showed that lensing and Hubble expansion can be modeled 
rigourously in a Modified Relativity. The degeneracy of MOND with Dark Matter can be broken if we examine 
the non-linear MONDian Poisson's equation in detail. Here we study the effect of tides for a binary stellar system 
or a baryonic satellite-host galaxy system. We show that the Roche lobe is more squashed than the Newtonian 
case due to the anisotropic dilation effect in deep-MOND. We prove analytically that the Roche lobe volume 
scales linearly with the "true" baryonic mass ratio in both Newtonian and deep-MOND regimes, insensitive to 
the modification to the inertia mass. Hence accurate Roche radii of satellites can break the degeneracy of MOND 
and dark matter theory. Globular clusters and dwarf galaxies of comparable luminosities and distances show a 
factor of ten scatter in limiting radii; this is difficult to explain in any "mass-tracing-light" universe. The results 
here are generalizable to the intermediate MOND regime for a wide class of gravity modification function ^{g) 
(Zhao and Tian, astro-ph/0511754). 
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1. Introduction 

The alternative gravity theory of Modified Newtonian 
Dynamics (MOND) (Milgrom 1983) has been doing very 
well in fitting kinematic data on galaxy scales, often bet- 
ter than the standard cold dark matter theory. Baryonic 
matter alone is sufficient to account for the gravity in 
such theory. The predictive power of this 20- year-old clas- 
sical theory with virtually no free parameters (Bekenstein 
& Milgrom 1984) is recently highlighted by the aston- 
ishingly good fits to contemporary kinematic data of a 
wide variety of high and low surface brightness spiral and 
elliptical galaxies; even the fine details of the ups and 
downs of velocity curves are elegantly reproduced without 
fine tuning of the baryonic model (Sanders & McGaugh 
2002, Milgrom & Sanders 2003). Originally it was pro- 
posed empirically (Milgrom 1983) that rotation curves of 
axisymmetric disk galaxies could be fit by an acceleration 
g = \g\ = V'^/r which is stronger than the Newtonian 
gravity GM/r'^ by a spatially varying factor 1/fi in the 
weak regime defined by g < oq ~ 1.2 x 10~^cmsec^^; e.g., 
/i(g/ao) = min(l,g/ao) < 1. This empirical MOND rela- 
tion can be elevated to a theory for an arbitrary baryon 
density distribution p(R), where a curl-free gravity field 
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g = — V$ is the gradient of a conservative potential $(R) 
and satisfies an equation (Bekenstein & Milgrom 1984) 
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Here we made the analogy with the Poisson's equation of 
a curl-free electric field E — —Wcj) generated by a cloud 
of static eletrons of density — e(R) in a medium with a 
spatially varying dielectric constant e{E) if we map the 
field g ^ E, the density p(R) ^ e(R) and the factor 
^ e{E). The above modified Poisson's equation is 
actually a reformat of the Lagrange equation ^ = —V • 
where the Lagrangian is given by 



- L= dK 
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g(R) = -V$(R), (2) 



where we integrate over the source energy density plus 

the field energy density; the factor fi{g) = i^!^^yi£j._ This 
classical field theory formulation guarantees conservations 
of energy and (angular) momentum of an isolated system. 

In the past this non-relativistic formulation of MOND 
has been criticized for being incomplete for modelling the 
bending of light (but see Qin, Wu & Zou 1995). This, too, 
has changed since its generalization into a respectable rel- 
ativistic theory (christened TeVeS by Bekenstein 2004), 
which includes Hubble expansion, and passes standard 
tests to check General Relativity (Skordis, Mota, Ferreira 
et al. 2005, Chiu, Ko & Tian 2005); GR is merely one 
limitting case of TeVeS. 
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Nevertheless, a main challenge of working on MOND 
is its essential subtle non-linearity and scale-dependency, 
which makes it unreliable to extrapolate Newtonian intu- 
itions. As a result, there are very few predictions of MOND 
in the literature in dynamical situations where the non- 
sphericity of the potential is essential. It is encouraging 
that the recent work of Ciotti & Binney (2004, see also 
Zhao 2005) shows surprisingly simple analytical scaling 
relations exist even for the highly non-linear and non- 
spherical two-body relaxation problem in MOND. Here we 
show a surprisingly simple scaling of tides or the Roche 
lobe of a binary system (on either stellar or galaxy scales) 
if it is in the non-linear deep-MOND regime. A subtle 
difference from a naive Newtonian extrapolation is also 
pointed out. We compare the predicted Roche lobe sizes 
with the observed limiting sizes of Milky Way satellites 
(globular clusters and dwarf galaxies) of 10^~^-^Lq. 

2. Roche lobe & binary potential in deep-MOND 

One way to reach the deep-MOND regime so that ^{g) — 
is to be at a distance R sufficiently far way from an iso- 
lated galaxy of total baryonic mass M so that the gravity 
g{R) <C ag. Here g{R) and the spherical galaxy potential 
^o{R) are approximately related to the Newtonian gravity 
GM/i?2 by 
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Consider introducing a low-mass satellite of mass m at 
a position {x, y, z) = (0, 0, Do) in the above galactic poten- 
tial ^o{R). Following Milgrom (1986), we consider the ef- 
fect of the self-gravity of a satellite mass inside an external 

galaxy field g{R) ~ \f^^^- The spatially slow- varying 
external field dominates the satellite gravity sufficiently 
far away from the satellite. So the MOND "dielectric in- 
dex" ^ ~ ii{g) ~ g{R)/ao ^ \f^^ varies very httlc in 
the vicinity of the satellite. To the first order in m the 
perturbation in potential is given by 
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where m' is the modified inertia of the satellite due to the 
external field go^ = giDo), and f is the effective distance 
from the centre of the satellite given by 
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g=g{Da) 



where 1 -I- A{g) is a shape factor; ^ clearly 1 + A{g) = 2 in 
the deep-MOND regime where /i ~ g/ao, and 1 + A{g) = 1 
for strong gravity. 

^ The "dielectric index" ^(■'f^) is more sensitive to pertur- 
bation along the external field g{Do)z direction than perpen- 
dicular because |g + dg\ = [(gOo + dg^)"^ + (dgx)^ + (dffy)^] ^ 
depends on the perturbation dg^ to first order, and dgx and 
dgy to second order. 
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Fig. 1. shows the re-scaled Roche lobes (contours of the 
effective potential) in the equatorial xz plane (lower half) 
and in the vertical yz plane (upper half) of a hypothetical 
isolated Earth-Sun binary with a mass ratio 3x 10~^ in the 
strong gravity regime (say with the separation Do = lAU, 
thin blue lines), and in the weak gravity regime (say with 
separation Do = O.lpc and GA/qD^^Oq ^ ~ 0.1, shaded 
areas) . The Earth is at origin and the Sun is at unit length 
to the right. The inner Lagrangian point is a saddle point 
between the Earth and the Sun, which is slightly further 
away from the Earth in the deep-MOND regime than in 
strong gravity regime. 

The above formulation allows us to approximate the 
potential of, e.g., the Milky Way galaxy with a satellite. 
Substitute in the expressions for /i, g{R), R and f, the 
combined potential is then given by 

$ = <S>o(R) + $1 = (GA/ao)5 lnv/z2 + + ^ (g) 

r 

which is an axisymmetric prolate potential with two cen- 
tres separated by distance Do along the z-axis, where the 
two terms represent the MONDian potential of the Milky 
Way and perturbation due to the satellite; here 
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We also note that the density p ~0 along the z-axis for the 
gravity field g = — V3> of the above combined potential; 
this can be verified by evaluating the modified Poisson's 
equation p — —V • 

Let the low-mass satellite with m/M ^ 1 rotate 
around the galaxy centre (fixed) with an angular veloc- 
ity rioYj then particles in the corotating frame conserve 
the Jacobi energy with an effective (triaxial) potential 



^e{x,y,z) = ^(x,y,z) 



'-^1 - (8) 



The inner or outer Lagrangian points is then calculated 
from the saddle point of the effective potential where 



5$e(0,0,z) 



dz 
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which defines the Lagrange radius r^. Taylor-expand the 
above to first order in tl/Dq (the 0th term cancels due to 
eq. |S1 and eq. [Tj) , we have 



M 

m 



i + C 



1 



So inside the Lagrange radius the average density of the 
satellite equals twice the average density inside the orbit 
of the satellite in the weak gravity regime. Note that the 
masses here m and M are true baryonic masses of the 
binary, not the modified inertia masses^. The scaling that 
r (X m^/^ is also confirmed by the numerical simulations 
of Brada & Milgrom (2000). 

The shape of the Roche lobe is defined by the con- 
tour of the effective potential (eq. |HJl passing through 
the Lagrange point. Finding the roots analytically yields 
(Zhao & Tian 2005) 



1, iVTo+V2)i- {VTo-V2)i 
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which are intersections with the long z-axis, intermedi- 
ate X-axis and rotation y-axis respectively. Of the three 
radii, the Intermediate Roche (IR) radius compares most 
directly with observed size in the sky plane for a distant 
satellite, and is given by 



Do 



V2 
3 



— ^ J_ _ = 0.374 



(12) 



Projection effects make the observed radius in between 
the short semi-axis and the long semi-axis of the Roche 
lobe. Hence the intermediate axis (nrt instead of r^) is 
the best compromise among the three to approximates 
the observed size. 

3. Observed instantaneous Roche lobe 

If MOND is correct the Roche lobe would act as Nature's 
balance to weigh the relative baryonic content of a sec- 
ondary vs. a primary star, or a satellite vs. its host 
galaxy. Interestingly the Roche lobe satisfies the same 
scaling relation ^ (f-)"^^^ = est, but the est = 0.462 
in strong gravity regime (Binney & Tremaine 1987) while 
est = 0.374 in deep-MOND. E.g., in a gedanke experiment 
where we take the solar system out of the Galaxy, and in- 
crease the Earth-Sun distance from lAU to O.lpc (the sep- 
aration of the widest known binary stars) so that near the 
inner Lagrangian point of the system the gravity drops 
from the strong regime to the weak regime. Fixing the 
Earth-Sun mass ratio m/M = 3 x 10~^, the rescaled inter- 
mediate Roche lobe radius tir/Z should decreases slowly 
by a subtle amount from 0.462 x (to/M)!/^ = 0.0067 ra- 
dian (for strong gravity) to 0.374 x (to/M)!/^ = 0.0054 

^ Dimensional analysis only cannot tell whether the dimen- 
sionless scales like (■p-)"or^-^^ > where m' = ^ is the 
modified inertia (cf. eq.|31l. 



radian (for weak gravity); cf. eq. (|10|1 and H12|l and see 
Fig.Q] Likewise the aspect ratios of the Roche lobe evolves 
from 1:2/3: Q^/^ _ 31/3 ^ 1 ; 0.667 : 0.638 (Binney & 
Tremaine 1986) to about 1 : 0.471 : 0.456, and the volume 



(10) of Roche lobe evolves from ' 



4irD„ 



to - 



47r£)„ 



The 
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Roche lobe is more squashed in MOND than in Newtonian 
gravity (cf. Fig.l). 

Interestingly, the same rescaled Roche radius can be 
predicted if we substitute the Earth-Sun binary by a satel- 
lite (either a dwarf spheroidal or a globular cluster) of a 
typical luminosity ~ 3 x lO^L© orbiting a luminous host 
galaxy of ^ lO^^L© so that the baryonic mass ratio is 
about Earth-Sun mass ratio. The self-gravity around an 
extended object of mass distribution m(r) becomes weak 
compared to ag outside a radius 



'Gm(r) 




m = IOMq, 
m = IO^Mq, 
10"Mp 
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Consider globular clusters and dwarf galaxies of the Milky 
Way much further than lOkpc. The outer envelope (well 
outside lOpc) of these objects are generally in the mildly- 
weak to the deep-MOND regime. A satellite is generally on 
a non-circular orbit, nevertheless, an instantaneous Roche 
lobe radius can still be defined by the hr as if the satellite 
is orbiting on a circular orbit at its present orbital radius 
Do (approximately the distance from the Sun D) for an 
outer halo satellite. So we can rewrite ea. 1121 as 



D \ Lmw J 



A = 0.374^ - 0.374, (14) 



where in estimating the instantaneous Roche radius riR 
(cf. ea. ll2|l we have assumed satellites have identical mass- 
to-light ratio as the Milky Way. This suggests comparable 
sizes riR for distant satellites at similar present distances 
D or Do and comparable luminosity Lgat- 

The actual direct observable is the limiting angular 
size 0iim of a satellite seen from the Sun's perspective. If 
a satellite fills the MONDian Instantaneous Roche Lobe, 
we expect to observe an angular size — tir/D per- 
pendicular to the line of sight. From these observables we 
can construct an observable "filling factor" 



rmjD 



A 



(15) 



In the real world, the limiting radius should be somewhat 
smaller than the instantaneous Roche radius, i.e, i^iim < 1 
because (i) the satellite likely remembers the truncation 
set by the stronger tide at some smaller orbital radius be- 
tween last pericentre and present distance D, depending 
on whether a satellite can relax quickly in one orbital time 
(Bellazzini 2004); (ii) this truncation might be beyond the 
limit of reliable observation riim simply because we run out 
of bright stars (Grebel, Odcnkirchen & Harbeck 2000). 
Many satellites (cf. Dinescu, Keeny, Majewski et al. 2004 
and references therein) are on nearly circular orbits (Ursa 
Minor and Ganis Major dwarves) or are presently within a 



4 



HongSheng Zhao: Roche Lobe Sizes in Deep-MOND Gravity 



factor of two to their pericentres (Pal3, Pall3, Fornax, Sgr 
and LMC). Theoretically a particle is typically found at 
the geometric mean ^Jryr^ of its pericentre and apocentre 
in a logarithmic potential. So D ^ Do ^ (1 — 2)rp typically 
even for a very radial orbit with Tp : ~ 1 : 4, i.e., we 
expect a mild scatter of i^nm between 0.5 — 1, allowing for 
radial orbits. The filling factor should be nearly unity for 
globular clusters since many galactic and some extragalac- 
tic globulars (Harris, Harris, Holland et al. 2002) appar- 
ently fill their Roche lobes judging from extra-tidal stars 
in power-law profiles revealed by deep observations wher- 
ever available (e.g., Leon, Meylan & Combes 1996). 

These expectations, however, are not borne out by 
Fig. [3 which shows surprisingly large scatter of the ob- 
served Fiim 0.1 — 5 for distant Milky Way satellites with 
comparable luminosities (10^ — 10^'^Lq). It is also difficult 
to understand why our expectation F\i-^ < 1 (cf. eq. I15|l 
is contradicted most strikingly by systems of larger tidal 
radii (larger symbols). Surely observations carry errors. 
The distance factor A is insensitive to the typical 10% 
distance error. Satellites often change profiles at 6'iim, so 
^lim is well-defined with very little error. Finally many 
satellites are in mild MOND regime with a Roche lobe 
size more rigorously given by (Zhao & Tian 2005) 



10.0 



Intermediate Roche Size 



1 + A TO 



Orbital Radius 3\/TTA V3 + A Af / ' ^^^-^ 

which is a smooth interpolation of eq. H12|l in deep- 
MOND and the familiar Newtonian prediction (Binney & 
Tremaine 1987) as the dilation factor (cf. eq.EJ 1 + A(^) 
changes from 2 to 1. For any /i- function of the Milky Way 
(cf. Famaey & Binney 2005), this correction is at most 
only 20% (the factor 0.374 in eq. US changed to 0.462). 

In short, it is likely challenging for any theory of struc- 
ture formation of a baryonic MOND universe to explain 
the puzzling large scatter in the rescaled Roche radius 
without fine-tuning of satellite orbits and mass-to-light 
ratios. It is less challenging for dark matter theories; 10^- 
star satellite objects could form either inside or outside a 
small dark halo. The scatter of satellite sizes echos with 
similar scatter of Einstein ring sizes around high-redshift 
lens galaxies (see lensing models of Zhao, Bacon, Taylor 
& Horne 2005), highlighting possible difficulties of mass- 
trace-light models. 
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